Mon. Not. R. Astron. Soc. 000,[T]{T6](2009) Printed 24 January 2010 (MN LOgX style file v2.2) 



Pressure Support vs. Thermal Broadening in the Lyman-a Forest I: 
Effects of the Equation of State on Longitudinal Structure 

Molly S. Peeples 1 *, David H. Weinberg 1 , Romeel Dave 2 , Mark A. Fardal 3 , Neal Katz 3 

1 Department of Astronomy and the Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210 

2 University of Arizona, Steward Obsen'atory, Tucson, AZ 85721 

3 'Department of Astronomy, University of Massachusetts, Amherst, MA 01003 



24 January 2010 



ABSTRACT 

In the low density intergalactic medium (IGM) that gives rise to the Lyman-a forest, gas tem- 
perature and density are tightly correlated. The velocity scale of thermal broadening and the 
Hubble flow across the gas Jeans scale are of similar magnitude (HXj ~ 0th)- To separate 
the effects of gas pressure support and thermal broadening on the Lya forest, we compare 
spectra extracted from two smoothed particle hydrodynamics (SPH) simulations evolved with 
different photoionization heating rates (and thus different Jeans scales) and from the pressure- 
less dark matter distribution, imposing different temperature-density relations on the evolved 
particle distributions. The dark matter spectra are similar but not identical to those created 
from the full gas distributions, showing that thermal broadening sets the longitudinal (line- 
of-sight) scale of the Lya forest. The turnover scales in the flux power spectrum and flux 
autocorrelation function are determined mainly by thermal broadening rather than pressure. 
However, the insensitivity to pressure arises partly from a cancellation effect with a sloped 
temperature-density relation (T oc p° - 6 in our simulations): the high density peaks in the 
colder, lower pressure simulation are less smoothed by pressure support than in the hotter 
simulation, and it is this higher density gas that experiences the strongest thermal broadening. 
Changes in thermal broadening and pressure support have comparably important effects on 
the flux probability distribution (PDF), which responds directly to the gas overdensity dis- 
tribution rather than the scale on which it is smooth. Tests on a lower resolution simulation 
(2 x 144 3 vs. 2 x 288 3 particles in a Mpc comoving box) show that our statistical 

results are converged even at this lower resolution. While thermal broadening generally dom- 
inates the longitudinal structure in the Lya forest, we show in Paper II that pressure support 
determines the transverse coherence of the forest observed towards close quasar pairs. 

Key words: cosmology: miscellaneous — cosmology: theory — intergalactic medium — 
methods: numerical 



1 INTRODUCTION 

The Lyman-a forest, caused by Lya absorption of neutral hydro- 
gen atoms along the line of sight to some distant source (usu- 
ally a quasar), was ori ginally described in terms of disc rete in- 
tervening gas "clouds" dLvndsl \l91 ll : ISargent et al.lll980l) . analo- 
gous to clouds in the Galactic interstellar medium. In this picture, 
the velocity widths of the observed absorption structures were pri- 
marily a consequence of thermal motions of the absorbing atoms. 
In the mid-1990s, three-dimensional hydrodynamic simulations of 
cold dark matter (CDM) cosmological models achieved remark- 
able s uccess in reprod u cing the observed properties of the Lya 
forest dCen et all 1 19941: IZhang et al.lll995l: IHernquist et al.|[l99r3 : 
iMiralda-Escude et al.lll996l : iTheuns et al.lll998l ). In these simula- 



tions, most of the absorbing gas is at low density (p/p ~ 0.1- 
10) and naturally described as a continuous ly fluctuating medium 



rather than a series of discrete structures (IHernquist et al 



iBi & Davidser]! 19971 : iRauch et alj[l997l : ICroft et all 1 19971 . 1 19981) 



19961: 



The velocity width of individual features is set largely by the Hub- 
ble flow across them, reflecti ng the physical extent of the absorbing 
gas a long the line of sight dHernquist et al.|[l996l : IWeinberg et al.l 
1 1997b . The temperature of the intergalactic medium (IGM), there- 
fore, affects the structure of the Lya forest in two ways: by smooth- 
ing absorption along the line of sight through the thermal motions 
of atoms, and by smoothing the physical distribution of the gas in 
three dimensions through pressure support. In this study, we use 
two smoothed particle hydrodynamics (SPH) simulations with dif- 
ferent photoionization heating rates — and thus different IGM tem- 
peratures and amounts of gas pressure support — to disentangle the 
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relative effects of pressure support and thermal broadening in the 
Lya forest. 

Both hydrodynamica l simulations dKatz et al.l 1 19961 : 
iMiralda-Escude et al.l 1 19961: I Theuns et al. I Il998h and analytic 
arguments i Hui & Gnedinl 19971) suggest that low-density inter- 
galactic gas should have a power-law temperature-density relation, 
i.e., T = T (l + S) a , where 1 + 5 = p/p is the local gas 
overdensity. Reasonable assum ptions for heating and cooling rates 
yield T ~ 10 4 K and a ~ 0.6 faui & Gnedinlll997l : lTheuns et all 
Il998h . Observations, however, imply that the normalization To 
could be nearly twice as high, and the slope could be much 
shallo wer or even inver t ed dSchave et al ] | 19991: 1 Ricotti et al ] |200d : 
iMcDonald etal] 1200 ll : iBolton et al.l l2008h . Regardless of the 
parameter values, because the optical depth to Lya absorption is 
related to both the temperature and the density by a power law, 
the existence of the temperature-density relation (also called the 
"equation of state") implies a tight relation between the observed 
Lya absorption and the IGM gas density. Since the universe has 
~ 5 times as much mass in dark matter as in baryons, we can 
expect in general for intergalactic gas to trace the underlyin g 
dark matter on scales above the gas Jeans length dSchavdl200lh . 
The Lya forest therefore provides a powerful tool for tracing the 
dark matter power spectrum in the quasi-linear regime — modulo 
the effect s of peculiar vel o cities, thermal broadening, and gas 
pressure llCroft et alj 1 19981, 1 19991. 120021: IMcDonald et all 1200(1 
l2006l:IViel et alj2004IViel & Haehneltll2006l) . 

Pressure should be important on scales at or below the Jeans 
length, 
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where c s = ^/[5fcT]/[3ra] = <r t h\/5/3 is the speed of sound in 
an ideal gas expressed as a multiple of the 1-D thermal velocity g t h 
jMiralda-Escude et alj [19961 : ISchavell200ll : iDesiacques & NusseJ 
120051) . Here m is the average mass of the gas particle; in an ion- 
ized primordial mixture of hydrogen and helium, m — 0.59m p , 
where m p is the proton mass. The density p is the density of the 
gravitating medium, which is dominated by dark matter, so we take 
(l + z) :, (l + <5).In comoving coordinates, 
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where we have normalized the thermal broadening velocity <T t h to 
correspond to a fiducial temperature of 10 4 K. By defining a "Jeans 
velocity" as vj = HXj we can find the relative importance of the 
Jeans scale and <7 t h, 
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which is essentially redshift-independent for the redshifts rele- 
vant to the Lya forest. The Jeans length of equation (QJ di- 
vides stable from unstable modes in a static, homogeneous, self- 
gravitating medium. The IGM is expanding, inhomogeneous, non- 
self-gravitating (because dark matter dominates), and evolving in 
density and temperature on the same timescale that fluctuations 
grow. Even for linear perturbations in a baryonic universe, the "fil- 



tering scale" below which fluctuations growth is suppressed de- 
pends on the thermal history o f the gas rather than th e instantaneous 
temperature-density relation jGnedin & Huilfl998l) . We therefore 
expect equation Q} to describe the scale of gas pressure support 
only at an order-of-magnitude level. Equation <[3j shows that Jeans 
velocities and thermal velocities should be comparable at the over- 
densities of typical Lya forest features, but the calculation is not 
definitive enough to show whether one will dominate in practice. 

Hence, to predict the statistical properties of the Lya forest 
in a given cosmological model, one must calculate the predicted 
gas distribution. The most reliable way to do this uses full N- 
body plus hydrodynamic simulations, which include the gravity of 
dark matter and gas and the additional effects of pressure, adia- 
batic heating and cooling, heating by photoionization and shocks, 
and radiative cooling. However, large volume hydrodynamic sim- 
ulati ons with the necessary reso lution are comput ationally inten- 
sive. IWeinberg et all i ll 9971) and ICroft etal] i 19981) show that one 
can achieve reasonable accuracy in the Lya forest regime from 
pure dark matter simulations, applying the temperature-density re- 
lation to the evolved dark mat ter distribution to com pute spectra. 
Indeed, the log-normal model dBi & Davidsen|[l997l) . in which the 
dark matter distribution is comp uted from the linear d ensity field 
by an exponential transformation dColes & Jones|[l99ll) , provides a 
qualitatively accurate physical model of the Lya forest, sufficient 
for creating artificial spectra with reasonable statistical properties. 

A fast way to include the effects of gas pressure in an ap- 
proximate way is the hydrodynamic part i cle-mesh (HPM) metho d 
(lGnedin& Hull 19981 : iRicotti et alJbood : | Mciksin & Whitdl200lh . 
HPM assumes that all the gas follows the power-law IGM equation 
of state, and it uses a modification of the standard (fast) particle- 
mesh iV-body method to compute the sum of gravitational forces 
and pressure gradient forces given this equation of state. In a de- 
tailed comparison of fully hydrodynamic SPH simulations and the 
approximated H PM simulations, both evolved using GADGET-2, 
IViel et all d2004l) found that while the HPM approach does con- 
verge to the SPH results, for some Lya forest properties, such as 
the flux probability distribution or small-scale power spectrum, it 
can differ from the SPH calculation by as much as 50% at z ~ 2. 
A simpler alternative to HPM is to use a pure TV-body simula- 
tion but smooth the evolved dark matter distribution on the Jeans 
scale before extracting Lya forest spectra fZaldarriag a et alJboOll : 
IDesiacques & Nusse J2005h . Because the Jeans smoothing is three- 
dimen sional, it is not degenera te with line-of-sight thermal broad- 
ening. [Zaldarriaga et alj feOOll) found that, even though the ther- 
mal broadening dominates the pressure correction, the value of 
the Jeans length becomes a large source of uncertainty in cosmo- 
logical inferences from the Lya forest if one tries to estimate it 
from the data rather than predict it from theory. Because the Jeans 
length is expected to vary with the gas temperature and density, a 
Zel'dovich-like scheme can be used when smoot hing the dark mat - 
ter distribution in order to model these effects dViel et al]l2002h . 
Though we do not carry out a comprehensive comparison of these 
methods here, we do investigate the impact of pressure in detail, as 
well as compare full SPH simulations to results derived from the 
dark matter distribution alone. 

The primary purpose of this paper and its companion (Peeples 
et al. 2009, hereafter Paper II) is to disentangle the roles of pressure 
support and thermal broadening in the Lya forest by studying two 
SPH simulations with different thermal histories and temperature- 
density relations. In addition to examining the physics of the Lya 
forest, this study is motivated by observational evidence (discussed 
in § 13. lb that the temperature of the IGM at z ~ 2—4 is higher than 
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expected from simple photoionization models by a factor of 1.5- 
2. This evidence comes from analyses of data along single lines 
of sight. Because the higher pressure associated with hotter gas 
would smooth the IGM in three dimensions, using closely paired 
lines of sight to probe this coherence scale has been proposed as 
an alternative route for inferring the temperature-density relation 
(J. Hennawi, private communication, 2007). However, before we 
can understand the relative roles of temperature and pressure on 
the transverse structure of the Lyman-a forest, we must first under- 
stand their longitudinal effects along independent sightlines, which 
is the goal of this paper. While here we find that thermal broaden- 
ing dominates pressure support in setting the level of longitudinal 
structure in the Lya forest, in Paper II we show that the gas Jeans 
length dominates the level of coherence transverse to the line of 
sight. 

This paper is organized as follows. In §|2] we describe the 
SPH simulations we have evolved to investigate these effects. In 
§[3] we describe the physics of the Lya forest in these simulations 
and examine the impact of the thermal history on observable spec- 
tra; this section also serves to review the physical understanding of 
the high-redshift Lya forest that has emerged from simulations and 
associated analytic work since the mid-1990s. We then study these 
effects on several typical statistical measures for learning about the 
IGM from the Lya forest in §|4] with our conclusions in §[5] 



2 SIMULATIONS 

We analyze two SPH simulations with identical initial conditions, 
one with a fiducial photoionization heating rate and one with a 
heating rate from photoionization that is four times higher than 
the fiducial; hereafter, we refer to these simulations as the "fidu- 
cial" and "H4" simulations, respectively. In the terminology of 
iKatz et al. I< fl996h . we compute the photoionization rates T and pho- 
toi onization heating r ates e for the fiducial simulation assuming 
the lHaardt&Madaul d200lh quasar + galaxy photoionizing back- 
ground, and for the H4 simulation we increase eh i, tHe i, and £He n 
by a factor of four. In principle, these heating rates could arise from 
a much harder UV background spectrum that yields more resid- 
ual energy per photo-electron. However, we do not propose any 
specific model for these heating rates — they are a computationally 
simple way to obtain IGM temperatures that are higher than those 
in the fiducial model and closer to those estimated from observa- 
tions (see 5 13, ll for further discussion). 

W e evolve these S PH simulations using the parallel GADGET- 
2 code dSpringellrlool) to trace the evolution of 288 3 dark matter 
and 288 3 gas particles in a 12.5/i _1 Mpc comoving cubic volume 
from z = 15 to z = 2. We also evolved another simulation using 
the fiducial heating rates and the same initial conditions but only 
2 x 144 3 particles to serve as a test for resolution convergence. We 
adopt a standard ACDM cosmological model with the parameters 
{QM,Q. A ,tt b ,h,a 8 ,n a ) = (0.25,0.75,0.044,0.7,0.8,0.95), all 
of which are in good agreement with the Wilkinson Microwav e 
Anisotropy Probe (WMAP) five-year results dHinshaw et al.l2009h ; 
our choice of o$ = 0.8, however, is somewhat lower than what 
is typically considered for Lya studies. These parameters give a 
mass per SPH particle of 1.426 x 1O 6 M0, which is much less than 
the typical Jeans mass of Mj = pA 3 7 ~ 7 x 10 9 Mq. The spline 
gravitational force softening has an equivalent Plummer length of 
0.875 h^kpc comoving (~ 1 /50 of the initial particle grid spac- 
ing). The SPH smoothing lengths are chosen to enclose 33 ± 2 
neighbors within the smoothing kernel. The simulation incorpo- 



rates standard heating and atomic cooling processes and the stan- 
dard GADGET-2 treatment of star formation and metal enrichment. 
These simulations do not incorporate galactic winds, but these 
should have very little effect on the Lya forest dKollmeier et al.l 
120061 ; iBertone & Whitel2006l ; iMarble et alj|2008l) . 

We extract spectra from the SPH gas distribution using 
TIPS^E as described by iHernguist et all d 19961) . Following com- 
mon practice, we rescale the intensity of the UV background so 
that the mean flux decrement of the extracted spectra matches ob- 
servations (Table [T), as the mean decrement itself is much bet- 
ter known than the background inten sity (see discussions by, e.g., 
ICroft et alj|2002l ; IMarble et al]|2008h . Our extracted spectra have 
1250 pixels across the 12.5/i -1 Mpc volume, making the pixel size 
w lkms -1 , which is well below the smoothing scale imposed by 
thermal broadening. 



3 PHYSICS OF THE LYMAN-a FOREST 
3.1 The Temperature-Density Relation 

In the absence of shock heating, the evolution of the temperature of 
the IGM is described by 
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as shown in detail by |Hui&GnedInldl997|) . Here, dS/dz and H(z) 
depend on the cosmology, the overdensity 1 + 8 = (p gas /pb), m 
is the mean particle mass, and dQ/dz is the net power per unit vol- 
ume owing to the ambient radiation field. The first two terms in 
equation (0 describe heating and cooling owing to adiabatic pro- 
cesses. After reionization, the change in temperature owing to the 
change in the ionization fraction (the third term in equation J4)) is 
effectively zero at all redshifts and relevant densities. 

Hydrogen reionization produces one energetic photoelectron 
per hydrogen atom, and it is expected to heat the IGM to a tempera- 
ture T ~ 2-5 x 10 4 K, depending on the s pectral shape of the ioniz- 
ing so urces and radiative transfer effects dMiralda-Escude & Reesl 
1 1994ft . Thereafter, adiabatic cooling reduces the overall tempera- 
ture, but denser regions remain hotter because they have higher 
neutral fractions and thus higher photoionization heating rates. At a 
given redshift, the simulated temperature-density (T-p) relation of 
the photoionized medium can be well approximated by a powerlaw, 



T{z) = T {z){l + 5) a ^\ 



(5) 



The slop e a(z) approache s 0.6 w ell after reionization. Evolution 
with the lHaardt & Madaul d200ll) UV background spectrum and 
fl b h 2 « 0.02 2 yields T sa 10 4 K at z = 2-4 for r eionization 
at_z> 7 (e.g. jHui & Gnedinlll997l ; lTheuns et alJll998l ; |Pave et al.l 
1 19991 ; Table [T). Figure Q] shows the distribution of SPH particles 
for our fiducial simulation in the temperature-overdensity plane at 
z = 3. At low overdensities (1+5 < 10), most of the gas falls along 
a tight locus, as expected from the above discussion. At higher den- 
sities (1 + 8 > 100), the gas has begun to cool to form galaxies. 
The higher temperature gas has been shock heated. (The apparent 
increase in temperature at 1 + 8 ~ 10 4 owes to the way GADGET-2 
treats the multiphase interstellar medium.) For many of our subse- 
quent analyses, we will isolate physical effects by imposing one of 
the three T-p relations, denoted by dashed lines in FigureQ] Specif- 
ically, we assign each gas (or dark matter) particle the temperature 



University of Washington version 



4 Peeples et al. 



implied by its overdensity and a given T-p relation before extract- 
ing spectra. The yellow line is an eyeball fit to the T-p relation in 
the fiducial simulation, while the pink line is the corresponding fit 
to the H4 relation; the H4 T-p normalization is ~ 2.3 times higher 
than the fiducial normalization. In both of these cases we set gas 
with 1 + S > 10 to a "shocked" temperature of T = 5 x 10 5 K, 
so that this high density gas will be entirely ionized and thus not 
contribute to the Lya forest. The mass- and volume-fraction of 
gas in these high density regions is relatively tiny; for all of the 
statistics presented here, the fiducial and H4 gas distributions us- 
ing the imposed fiducial and H4 T-p relations, respectively, yield 
nearly identical results to those obtained using the actual tempera- 
tures calculated by GADGET-2. Parameters for these fits are listed 
in Table [TJ For some comparisons, we impose a flat T-p relation, 
with all particles set to T — 2 x 10 4 K, as shown by the blue dashed 
line. 

Although most Lya forest features are broadened by Hub- 
ble flow, the narrowest features arise at velocity caustics and 
have widths set by thermal broadening. Analyses of observed 
line widt h distributions imply T sa 1.5-2 x 10 4 K at 2 < 
z ^ 4 dGnedin&Huilll998l ; iRicotti et alj |2000|; ISchaye et al.l 
l200d ; iMcDonald et alj l200ll ; we adopt iMcDonald et all s con- 
straints in Table [T] and Figure O. This is significantly hotter than 
the value T « 10 4 K expected for a lHaardt & Madaul feOOlh 
ionizin g background. iTheuns et al. I d200(jT and IZaldarriaga et^l] 
feOOlh find a similar result by fitting the small-scale cutoff of the 
one-di mensional flux power spectrum. More recently. iBolton et alj 
d2008l) have fit the flux probability distribution function (PDF) in 
high-resolution spectra inferring a similar, hig h Tp an d a shallow, 
possibly inverted (a < 0) slope. iLidz et all d2009i ) find To ~ 
2 x 10 4 K over the redshift range z = 2-4, using a wavelet analy- 
sis of 40 high-resolution spectra. In all cases, the observations are 
interpreted by comparing them to a suite of cosmological simula- 
tions, which incorporate many parameters in addition to the T-p 
relation itself. 

According to equation there are two basic ways to increase 
the gas temperature: either have hotter "initial conditions," or a 
higher heating rate dQ/dz. Even with the highest plausible reion- 
ization temperatures, it is difficult to reproduce the infe rred Tp at 
z = 3, given the observatio nal evid ence that -z r eion ^ 7 iFan et alj 
l2006l ; lHinshaw et alj|2009l ; see also lHui & Haimarj|2003l) . Energy 
injection from He II reionization over an exte nded period of time 
could also keep the IG M hot down to z ~ 3 l lBolton et a"l]|2008l ; 
iFurlanetto & Oh1l2008l) . though even this mechanism appears un - 
able to produce an inverted T-p relation dMcOuinn etaT]|2009l) . 
A harder photoionizing spectrum produces more energy input per 
photoelectron and thus a higher IGM temperature^ It is this solu- 
tion that we adopt for our H4 simulation, though the heating rates 
we adopt correspond to an implausibly hard spectrum. The struc- 
ture of the Lya forest should depend on the T-p relation but be 
fairly insensitive to the detailed mechanism that produces it, so we 
expect our conclusions about the impact of pressure support to ap- 
ply to a broad range of such mechanisms. However, the tension 
between the T-p relations predicted from theory and those inferred 
from observations remain puzzling, and it is not clear whether the 
resolution lies in modest changes (e.g., temperatures at the low end 
of observational estimates and a spectrum that is harder than con- 



2 Note that in photoionization equilibrium, higher intensity background 
radiation — with the same spectrum — affects only the photoionization and 
recombination rates, not the electron temperature. 
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Figure 1. Distribution of 1 % of the gas particles for the fiducial simulation 
in the temperature-density plane at z = 3. In subsequent analyses, we use 
either the simulation temperatures themselves or one of the three plotted 
T-p relations (dashed lines). 



ventionally assumed) or a physical process that has not yet been 
identified. Precisely because of this tension, it is important to un- 
derstand the physics and observational consequences of pressure 
support in the Lya forest, as we study in this paper and in Paper II. 

In Figure [2] the dashed lines show the approximate T-p 
relationships found in our two SPH simulations; the z — 3 
dashed lines are the same as the low-density yellow and pink 
lines in Figure Q] The temperature-density parameters for the ob- 
servations and our simulations at each redshift are given in Ta- 
ble Q] The fiducial simulation clearly disagrees with the obser- 
vations at z = 3 and 2.4, while the H4 simulation is in closer 
agreement but somewhat too hot, i.e., the two simulations bracket 
the central observation estimates. The points in Figure [2] are nu- 
merical integrations of equation <|4}, where we approximate the 
overdensity evolution using a modifie d Zel'dovich approximation 
dReisenegger & Miralda-Escudej [l995). H I reionization is mod- 
eled by initializing the temperature at a "reionization temperature" 
T r for all overdensities at a reionization redshift 2 re ion = 9.45; the 
temperature-density relation is fairly independent of T r by z = 4, 
assuming reionization occurs at z > 6. 

3.2 The Fluctuating Gunn-Peterson Approximation 

The H I optical depth r H i is proportional to the neutral hydrogen 
density, with 

(6) 



Tm = /Ao-ff (z)n m 



where fx is the oscillator strength of the Lya transition and Ao is 
the center of the Lya transition, 1216A dGunn & Peterson! [l9ol : 
lMiralda-Escudelll993l) . In the limit of high ionization fraction, the 
neutral hydrogen density, 
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Figure 2. A comparison of the temperature-density relation evolution from sim ulations {dashed l ines), calculations using equation (4} {points), and observations 
{shaded regions). The shaded region corresponds to the l-cr error region fromlMcDonald et al.1 f200ll) , with the orange lines indicating the best-fit relations at 
each redshift. (There are two such lines for z = 3 and 2.4 because|fylcDonald e t alj used simulations at two different redshifts to compare to these observations.) 
At z = 3, the dashed pink and black lines are the same as the dashed pink and gold lines, respectively, in Fig.[T] 

Table 1. Observed mean flux de crements (D) = (1 — e~ T ) are from lMcDonald et alj l2OO0l) and observed temperature-density relations (T = T)[l + 8] a ) 
are from lMcDonald et al.1 feOOll) . 
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agnostic of the IGM temperature To, but in practice it is not well 
enough known. Hence, we rely on structure in the Lya forest for 
IGM diagnostics and cosmological tests. 

Our results below rely on our full hydrodynamical simula- 
tions, but equation l[8) is useful to understand our results, and it 
can be a useful basis for simpler analytic or numerical treatments. 
It is often referr ed to as the "fluctuating Gunn-Peterson approxi- 
mation" (FGPA; I Weinberg et ail 1 19971: ICroft et al.|[l998n because 
it describes Lya absorptio n as a continuous phenomenon analo- 
i gous to the Gunn-Peterson {l965) effect, but arising in a fluctuating 
medium. As written, it ignores thermal broadening (a temperature- 
and therefore density-dependent convolution) and shock heating 
(i.e., gas not falling on the temperature-density relation). Also, 
while 1+5 = Pgas/pgas, the gas overdensity is often approximated as 
the dark matter overdensity Pdm/Pdm in TV-body simulations, per- 
haps smoothed by an effective Jeans length. 



3.3 The Lya forest in the fiducial case 

Figure[3]shows, on the left, a slice through our fiducial simulation 
at z = 3, with a depth of 125 hT 1 kpc. The small-scale filamentary 
structure of the high redshift universe is evident. While the dark 
matter and baryons have similar large-scale structure, the zoom-in 
panel on the right shows that the gas distribution is more diffuse. 
In particular, the densest filaments of dark matter lie within thicker 



is proportional to the gas density squared; Tuv is the photoioniza- 
tion rate owing to the ambient UV background. At the relevant tem- 
peratures (T <C 10 6 ), the recomb ination coef ficient a H n is propor- 
tion al to T~° MKatz et alj| 19961) . Following lilauch et alj Jl997h 
and lCroft et al.1 l ll998l) . wecan combine equations lO and 10 with 
a power-law T-p relation to obtain the optical depth 
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The last factor represents the impact of line-of-sight peculiar ve- 
locity gradients AV\ os /Ax, which change the density of atoms in 
frequency space relative to real space. 

Independently of the Lya forest, the photoionization rate Tuv 
is difficult to constrain; the other quantities entering the normal- 
ization pre-factor of equation ® a re also uncertain. We there- 
fore follow standa rd practice (e.g., iMiralda-Escude et al.l 1 19961 : 
iMarble et al1l2008l) and choose Tuv for each model so that it re- 
produces the observed mean flux decrement (D) = (1 — F) = 
(1 — exp(— T H i)} at the r edshift under investigat ion. Specifically, 
we adopt (D) values from lMcDonald et alj j200ol). which are listed 
in Table[TJ at the overlapping redshifts {z = 2.4 , 3) these ( D) ar e 
consistent with the more recent measurements of lKim et al.1 1 |2007|) . 
If Tuv were perfectly known, then (D) could itself be used as a di- 



6 Peeples et al. 




matter particles (black), and star particles (blue) all shown. The five green lines (left) denote sightlines referred to in subsequent figures. 



filaments of gas. This difference reflects the impact of gas pressure 
support. 



Figure [4] shows how density, temperature, velocity, and ther- 
mal broadening combine to produce Lya forest spectra, along the 
sightlines marked 1-4 in Figure [3] In each panel, the bottommost 
plot shows the gas temperature in black and the neutral hydrogen 
fraction in blue. Because n m /n H oc n H T~ ' 7 and T oc n^' 6 for 
gas on the T-p relation, the neutral fraction and temperature are 
positively correlated at most temperatures. However, high temper- 
ature gas has usually been shock heated off the T-p relation, and 
the recombination c oefficient q H ii fa lls more steeply than T~ 0,7 
at high temperature dKatz et aT]|l996f) , causing the neutral fraction 
to decrease dramatically (see, e.g., the feature near 100kms~ in 
sightline #2). In the middle graphs, we plot the neutral hydrogen 
number density [cm"' 1 ] in purple, using the scale on the left-hand 
axis. The gas and dark matter overdensities are plotted in black and 
cyan, respectively, using the scale on the right-hand axis. In gen- 
eral, the dark matter and gas overdensities are similar, but the dark 
matter has sharper, higher overdensity peaks as was seen visually 
in Figure [3] As expected from the T-p relation, the gas overden- 
sity and the gas temperature follow one another except at very high 
gas overdensity. The highest gas overdensities correspond to con- 
densed halos — i.e., galaxies — and thus the gas has cooled to lower 
temperatures in these regions. The neutral hydrogen density shows 
more variation than the gas density because it is proportional nf, 
(actually p 16 once temperature effects are included). The topmost 
plot has the transmitted flux F = exp(— r ra ) in black and the trans- 
mitted flux that would be observed in the absence of thermal broad- 
ening in grey. The lines between the middle plot and the top plot 
show the effects of peculiar velocities when converting from neu- 
tral hydrogen density in physical space to an observable flux in 



velocity space0 Features contracting along the line of sight cause 
physically distinct gas regions to converge to the same region of the 
observed spectrum. However, if the connecting lines do not cross, 
then the region still has net expansion. Both the thermally broad- 
ened and non-thermally broadened spectra have the same character- 
istic broad features, implying that residual Hubble flow dominates 
the velocity width of these features. However, thermal broadening 
smooths the small scale roughness. At velocity caustics (converg- 
ing lines in Figure @, peculiar velocities cancel the Hubble flow, 
and the features do become narrower when one removes the ther- 
mal broadening. 

3.4 Impacts of temperature and pressure on gas evolution 

The obvious consequence of having a higher photoionization heat- 
ing rate is that gas temperatures in the H4 simulation are higher 
than in the fiducial simulation, as shown at z = 3 in the top pan- 
els of Figure [5] A more subtle effect, is that the larger Jeans length 
of the H4 simulation smooths the gas distribution, as shown in the 
density-coded bottom panels of Figure[5] There is a relative paucity 
of very dense clumps in the H4 simulation; this is especially notice- 
able in the lower density filaments. 

Figure [6] plots the distribution of gas overdensity in the fidu- 
cial and H4 simulations. The visual differences in the right panel 
of Figure [5] manifest themselves as a higher frequency of particles 
with Pgas/pgas ~ 10-100 in the fiducial simulation and a higher 
frequency of p g as/pgas ~ 1 particles in the H4 simulation. This 
difference leads to a cancellation in many of the statistical compar- 
isons in §|4j though the H4 gas is hotter (and thus has more ther- 
mal broadening), its higher pressure implies that there there is rela- 
tively less high-density and therefore relatively higher-temperature 

3 The bar-like appearance of the lines showing the effects of peculiar ve- 
locity also lends these complicated plots the name of "zoo plots." 
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Figure 4. Physical quantities and Lya forest spectra along four of the sightlines labelled in Figure[3] velocity v = corresponds to the top of Figure[3] Shown 
in real space are the temperature (black, bottom), the neutral hydrogen fraction (blue, bottom), the H I number density (purple, middle), the gas overdensity 
(black, middle), and the dark matter overdensity (cyan, middle). The bars separating the top and middle panels show the effects of peculiar velocities when 
transitioning from real space to the observed flux transmission with (top, black) and without (top, grey) thermal broadening. See § l3.3l for more details. 
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Figure 5. Temperature (top) and density (bottom) evolution for gas particles 
in a 2.5 X 2.5 X 0.05 /i -1 Mpc comoving slice at z = 3 in the fiducial (left) 
and H4 (right) simulations. The temperature scale runs logarithmically from 
from logT = 3.5 (black/dark blue) to logT = 5.5 (yellow/white). The 
density scale runs from log(p sas /p gas ) = log(l + 5) = —1.3 (black) to 
log(l + 8) = 1.3 (green). 



gas. For comparison, we also show in Figure [6] the density distri- 
bution of the dark matter, with densities computed using the SPH 
smoothing kernel. Because the dark matter is pressureless, it typi- 
cally reaches higher overdensities, though it does not achieve the 
highest overdensities seen in the gas distributions because these 
arise from dissipation, i.e. cooling. Extrapolating from the density 
distributions in our low-resolution simulation, the dark matter den- 
sity distribution would probably be somewhat less skewed if we 
increased the mass resolution of the simulation, but at low overden- 
sity the gas distributions are essentially converged since we resolve 
the Jeans mass. 



4 EFFECTS OF PRESSURE AND THERMAL 
BROADENING ON THE Lya FOREST 

To isolate the effects of gas pressure and thermal broadening in 
our subsequent analyses, we examine both the original gas parti- 
cle distributions and distributions with one of the three imposed 
T-p relations illustrated in Figure [TJ For an imposed T-p relation, 
we replace each gas particle's temperature by the temperature that 
corresponds to its overdensity. We use the same procedure to cre- 
ate Lya forest spectra from the dark matter distribution ("gasify- 
ing" the dark matter). The fiducial T-p relation (yellow in Fig- 
ure [TJ matches that found for the IGM in the fiducial simulation; 
above 1 + S — 10, the gas is "shock heated" to a temperature of 

5 x 10 5 K. Likewise, the pink lines in Figure [TJ show the T-p re- 
lation used to mimic the H4 simulation. Figure [6] shows that gas 
distributions with different pressure also sample the T-p relation 
differently, and if the T-p relation is sloped they will therefore ex- 
perience different thermal broadening For example, using the fidu- 
cial T-p relation, the H4 simulation would have less high density 
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Figure 6. Distributions of gas overdensities at z = 3 in the fiducial (black) 
and H4 (pink) simulations, as well as the dissipationless dark matter (cyan) 
in the fiducial simulation. 



gas with large thermal broadening, and the dark matter distribution 
would have more. We, therefore, consider an additional, flat T-p 
relation, with T = 2 x 10 4 K at all densities, so that we can exam- 
ine the effects of pressure in the presence of pressure-independent 
thermal-broadening. 

We now turn to the effects of pressure support and the T-p 
relation on the Lya forest spectra (§ 14. It and on flux statistics (the 
flux power spectrum, § 14. 21 the autocorrelation function, § 14.31 and 
the probability distribution function, § 14.41 For all these analyses 
we use the same sightlines in each simulation to ensure the effects 
of sample are variance the same. 

4.1 Spectra 

Figure [7] shows how pressure and the T-p relation affect the Lya 
spectra along the five sightlines labelled in Figure [3] with sight- 
line #1 corresponding to the topmost spectrum in each panel; as 
in Figure [3] v = corresponds to the top of Figure [3] Any dif- 
ferences between the models are usually most noticeable in unsat- 
urated lines. The top-left panel shows spectra computed directly 
from the simulated gas distributions, with differences in pressure 
support, differences in thermal broadening because of the different 
T-p normalization, and differences in thermal broadening because 
of the sampling of the T-p relation. Light grey lines show spectra 
from our lower resolution simulation. The impact of resolution is 
generally very small, but there are some slight differences. 

There are two notable differences between the fiducial and H4 
spectra exemplified by the v ~ 450 km s -1 feature in sightline #5 
and the v ~ 750 km s -1 feature in sightline #4. First, features in 
the H4 spectra are broader than in the fiducial spectra, which could 
arise because of greater thermal broadening and/or because of the 
larger Jeans scale (and thus the larger HXj). Second, features in 
the H4 spectra are not as deep as in the fiducial case because the 
hotter gas does not reach as high overdensity and/or because higher 
thermal broadening smears out inherently sharp features. 

The top right panel isolates the impact of thermal broadening, 
applying the fiducial and H4 T-p relations to the gas distribution 
of the fiducial simulation. The lower-left panel isolates the impact 
of pressure, applying the flat T-p relation to the fiducial and H4 
gas distributions and to the dark matter distribution from the fidu- 
cial simulation. In some cases, such as the v ~ 150 km s _1 feature 
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Figure 7. Comparison of spectra along the sightlines labelled in Figure [5] isolating different physical effects. Top left: full effects of thermal broadening, 
thermal history and resolution are shown. Top right: effects of thermal broadening are isolated. Bottom left: effects of pressure support are isolated. Bottom 
right: effects of pressure support and different samplings of the underlying overdensity distribution are isolated. See § 14. 1 I for details. 
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in sightline #1, the v ~ 1025 km s -1 feature in sightline #3, and 
the v ~ 250 km s -1 feature in sightline #4, differences in the full 
simulation spectra (upper left) are largely erased in the upper right 
panel but remain similar in the lower left panel, which shows that 
they mostly arise from different pressure effects in the two simula- 
tions rather than from differences in the thermal broadening. There 
are fewer cases of the reverse, where differences present in the full 
spectra remain in the upper right panels but disappear in the lower 
left, though some are visible in sightline #5. The differences in 
the upper-right (thermal broadening isolated) are always system- 
atic, with lower thermal broadening yielding deeper and slightly 
narrower features. The differences between the H4 and fiducial gas 
distributions in the lower-left (pressure isolated) are more random, 
and we will see below that their statistical signature is weaker. 

The pressureless dark matter distribution does lead to spectra 
with more small scale structure, as shown by the thin dark lines in 
the lower left. However, despite the clear differences between the 
gas and dark matter distributions evident in Figures[3][4] and|6] the 
differences in the spectra created from these distributions are small. 
In the absence of thermal broadening, dark matter spectra are much 
more jagged than their gas counterparts shown in Figure[4] but real- 
istic thermal broadening masks these differences to a large extent. 
The velocity widths of most features are not sensitive to the level 
of thermal broadening (as seen in the upper right) but neither are 
they set by the Jeans scale, or else the dark matter and gas spectra 
would have greater differences. Instead, thermal broadening erases 
the finest scale structures in the density field so that typical features 
correspond to coherent, moderate overdensity structures expanding 
with the Hubble flow. 

Adopting the fiducial T-p relation in place of the flat T-p re- 
lation (bottom right) makes only a small difference to the usual ap- 
pearance of the spectra (top left). Features that are flat-bottomed 
in the dark matter spectra with the flat T-p relation often be- 
come less saturated with the fiducial T-p relation, such as the 
v ~ 775 km s -1 feature in sightline #3 and the v ~ 1000 km s~ 
feature in sightline #5. In these regions, the dark matter overdensity 
is high, so with T oc (1 + <5) 0,6 they have higher thermal broad- 
ening, which spreads the feature in velocity space and reduces its 
saturation. In some cases, such as the v ~ 725 km s -1 feature of 
sightline #1, the wings of the dark matter feature become notice- 
ably broader than those of the gas features because of the higher 
temperatures at higher overdensities. 

4.2 The 1-D Flux Power Spectrum 

The Lyman-Q flux power spectrum is a powerful tool for prob- 
ing the dark matter mass power spectrum on the smallest scales. 
Because of the high redshift, the continuous sampling of the line- 
of-sight density field, and the moderate overdensity of absorbing 
structures, the Lya forest provides a more direct link to the lin- 
ear theory power spectrum than other small-scale tracers. (Pri- 
mary cosmic microwave background anisotropies are damped on 
these scales.) However, to infer information about the structure 
of the underlying dark matter, one must understand the thermal 
structure of the IGM, as the amplitude and the shape of the flux 
power spectrum are connected to the gas temperature-density re- 
lation via equation l[8}. The bias of the flux power spectrum, i.e., 
b 2 (k) = Pf(fc) / Piin(k), is also influenced by non-linear gravita- 
tional evolution, thermal br oadening, pressure support, peculiar ve- 
locities, and shock heating dCroft et al'l 19981,1 1999l . l2002MViel etakl 
l2004l2008l ; lMcDonald et all 20051) . 

Figure [8] shows the one-dimensional line-of-sight flux power 
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Figure 8. The 1-D flux power spectrum at z = 4.0 (top), z = 3.0 (middle), 
and z = 2.4 (bottom), with the fiducial simulation in black, the H4 sim- 
ulation in pink, and the lower-resolution simulation in grey. For reference, 
we also show the flux power spectrum in the absence of thermal broadening 
for the fiducial gas (green) and the H4 gas (purple ), as well as t he observed 
flux po wer spectrum from McDonald et al. (200]], squares) and lCroft et alj 
l2002t triangles). 
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spectrum, Pio(k) of (F — (F)), where F = exp(— r Lya ), based on 
600 randomly selected sightlines through the fiducial, H4, and low 
resolution simulations. The flux power spectrum Pm(k) is a mea- 
sure of the variance of the flux, up, on diff erent scales; spec i fically , 
we adopt the normalization convention of iMcDonald et all d2000l) . 
where 
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Also plotted in Figure [8] are the flux power spectra for the same 
lines of sight in both the fiducial and H4 simulations in the ab- 
sence of thermal broadening. For reference, we also plot the ob- 
served flux power s pectrum measured from high-resolution spectra 
given in Table 4 of IMcDonald et al.l d2000L squares) and Table 7 
(the "B" sample at z = 2.4 and the "D" sample at z = 3) of 
ICroft e t aL (2002, triang les). We have ren ormalized the lCroft et al.l 
d2002h measurements bv lMcDonald et alT s {F} 2 to match our nor- 
malization convention. 

The cutoff in the flux power spectrum at large k is a conse- 
quence of thermal broadening, as is obvious from comparing the 
power spectra with and without thermal broadening. With thermal 
broadening included, the low resolution (144 3 gas particles) and 
fiducial (288 3 ) simulations produce very similar power spectra, in- 
dicating that even the lower resolution simulation is well converged 
for this statistic. At high k, the power spectra of the H4 simulation 
are offset by roughly a factor of 1.6 in k from the fiducial simula- 
tion, roughly consistent with the difference of 2-2.5 in To (since 
1/2 

thermal velocities scale as T ). The cutoff scale in the fiducial 
simulation agrees better with the observa tional data, even t hough 
its temperatures are low compared to the IMcDonald et alj fcOOll) 
estimates (see Fig. [2}- 

Our cutoff scale also differs from that of llidz et al. I J2009l, 
Figure 5); the flux power spectrum of our fiducial simulation is 
more similar (on small and large scales) to that of their hotter, 
To = 2 x 10 4 K and a = 0.3, simulation than to their colder, 
T) = 1 x 10 4 K and a = 0.6, simulation, which has a T-p rela- 
tion similar to our fiducial simulation. We see no obvious reason 
for this discrepancy: the resolution test in Figure [8] shows good 
convergence, we do not expect the cutoff scale to be sensitive to 
box size, we have checked that high as density peaks produce the 
expected thermally broadened line profiles in our extracted spectra, 
and we find good agreement between the flux power spectrum from 
our 2 = 3 spectra measured by our code (written by R. Croft) and 
an independent code written by P. McDonaldR When spectra are 
extracted from our simulation using iLidz et all s code, the calcu- 
lated P(k) is consistent with the one we measure the difference 
is therefore present in the simulated gas distributions themselves. 
For now, we draw no strong conclusions from the comparison to 
data in Figure(8] and focus instead on the relative roles of pressure 
support and thermal broadening. 

To separate the effects of gas pressure and thermal broaden- 
ing, Figure [9] shows z = 3 flux power spectra from 200 randomly 
selected sightlines through the fiducial gas, H4 gas, and fiducial 
dark matter density fields (indicated by the line type), each with 
three imposed temperature-density relations: fiducial, H4, and flat 
(indicated by the line color). Comparing lines of the same color 
in Figure [9] shows the effects of pressure support, with the same 
T-p relation applied to different density distributions. Comparing 



4 We thank Patrick McDonald for carrying out this test for us. 

5 We thank Adam Lidz for carrying out this test for us. 
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Figure 9. The 1-D flux power spectrum at z = 3 for fiducial gas (solid 
lines), H4 gas (dashed lines), and fiducial dark matter (dotted lines) with 
fiducial (brown), H4 (red), and flat (blue) imposed temperature-density re- 
lations. 



lines of the same type but different colors isolates the effect of the 
T-p relation for the same underlying density distribution. The lines 
clearly separate into three groups based on color, showing that ther- 
mal broadening dominates over pressure support in determining 
the scale of the power spectrum cutoff. There is some difference 
among the three density distributions when we impose a constant 
IGM temperature T = 2 x 10 4 K (blue lines), which shows that 
some of the similarity for the other T-p relations may reflect the 
cancellation discussed in § 13.41 where the distribution with a larger 
Jeans length has less high density gas to experience high thermal 
br oadening. Howeve r, on the whole our results confirm the finding 
of lMcDonaldl d2003l) that it is the instantaneous T-p relation (at the 
epoch of observation) rather than the detailed gas thermal history 
(and thus pressure history) that determines the Pf (k) cutoff. 



4.3 Flux Decrement Autocorrelation Functions 

The flux decrement autocorrelation function, 

(D(v)D{v + Av)} 



£auto(A«) 



(Dr- 



ew 



is a commonly used tool for characterizing the Lya forest. While 
technically the one-dimensional flux power spectrum and flux au- 
tocorrelation functions codify the same information (Pf (k) is just 
the one-dimensional Fourier transform of £ a uto), the flux autocorre- 
lation function is often easier to describe both observationally and 
from simulations because its definition does not depend on infor- 
mation from all scales. Furthermore, as cross-correlation functions 
are more typically used to describe information in closely paired 
lines of sight than cross-power spectra, it is important for us to 
understand the effects of pressure and temperature on the autocor- 
relation function before examining paired sightlines (in Paper II). 

Figure \10\ plots normalized flux decrement autocorrelation 
functions at z = 2.4, 3.0, and 4.0. We normalize £ au to = 1 at 
the smallest scales to highlight differences in the turnover scale of 
the correlation function rather than differences in normalization. 
The fiducial and low-resolution simulations have nearly identical 
Cauto, while the H4 simulation has a noticeably larger coherence 
scale. However, if we impose the fiducial T-p relation on the H4 
gas, then £ au to is nearly identical to that of the fiducial simulation. 
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Figure 10. Left: Normalized flux decrement autocorrelation functions, gauto = (D{v)D{v + Av))/{D) 2 , at z = 2.4, 3.0, and 4.0. The black curve with 
eiTor bars shows results for the fiducial simulation, while the pink and grey lines show results for the H4 and low-resolution simulations, respectively. Red and 
mustard lines show, respectively, the fiducial gas with the H4 T-p relation imposed and the H4 gas with the fiducial T-p. The green line shows the fiducial gas 
with no thermal broadening. Right: Normalized £auto, for the fiducial gas (light blue), H4 gas (bright blue), and fiducial dark matter (dark blue), all with an 
imposed temperature of T = 2 X f 4 K. 



Conversely, imposing the H4 T-p on the fiducial simulation yields 
nearly the same £ au to as the H4 simulation[f| 

While these results suggest that the effects of pressure are 
small compared to those of thermal broadening, the right-hand pan- 
els of Figure[l0]show that £ au to is different for the H4 gas, fiducial 
gas, and fiducial dark matter if we impose a flat T = 2 x 10 4 K on 
all particles. In this case, the higher pressure distribution exhibits a 
larger coherence scale, especially at high redshifts. Thus, the sim- 
ilarity of the fiducial and H4 distributions with the same imposed 
T-p, seen in the left-hand panels, arises partly from the counter- 
balancing effects of To and the changes in the distribution of gas 
densities discussed in § !3.4l (see Figure[6]l. 



6 If we impose the fiducial (H4) T-p relation on the fiducial (H4) sim- 
ulation, instead of using the simulation temperatures themselves, then the 
change in § auto is negligible. 



4.4 Flux Decrement Probability Distributions 

The flux decrement probability distribution function (PDF) is 
a potentially powerful tool for probing the IGM, both by it- 
self and as a means of adding constraints to other statistical 
measures such a s the flux power spectrum dRauch et al. 19971; 
Weinber gl 1 19991; iGaztanaga & Croftl 1 19991: iNusser & Haehneltl 



2000; iDesiacques & Nu sser 2005; iDesiacques et al.|[2007l) . In par- 
ticular, Bolton et all d2008l) " have recent l y used the Lya forest flux 
distribution as measured by iKim et ail 020071) to suggest that the 
temperature-density relation at z ~ 3 is inverted (a < in equa- 
tionO, with low-density gas at higher temperature than higher den- 
sity gas. While the differences owing to changes in the T-p relation 
are not as extreme as those in Pp(fc) or £ au to, the flux decrement 
PDF can potentially be measured to much higher accuracy using the 
same number of sightlines. On the other hand, the PDF is sensitive 
to continuum fitting, while on small scales Pf (k) is not (on large 
scales continuum fitting can be problematic even for Pp(k); e.g., 
I Kim et al. I l2004h . Furthermore, high-resolution spectra are neede d 
to accurately measure the PDF dViel et al.l2004l ; |Tvtler et al.l2004l) . 
Since high-resolution spectra will always exist in smaller numbers 
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Figure 11. Flux decrement probability distributions in linear bins of D at 
z = 2.4, 3, and 4 for the fiducial and H4 simulations, with (black and 
pink) inherent temperatures and imposed H4 (red) and fiducial (mustard) 
T-p relations, respectively. The points at z = 2.4 and 3 are included to 
show the p(D) in the most transparent (left) and opaque (right) bins. 



than low-resolution spectra, it should be noted that the PDF is also 
otent ially sensitive to sample variance, the explanation lKim et al.l 
20071) suggest for the ~ 30% difference between their z ~ 3 PDF 

and the one measured bv lMcDonald et al.1 fcOOOh on a very similar 

data set. 

Figure QT] shows the linear flux decrement PDF, i.e., p(D), 
where p(D)AD is the number of pixels in a bin of width AD di- 
vided by the total number of pixels. In Figure [T2] we show for a 
larger set of models the logarithmic PDF, p(logD), which allows 
better visual discrimination in the range 0.1 < D < 0.5. The shape 
of the PDF changes radically with redshift as the mean flux decre- 
ment changes; physically, this change is driven by the mean density 
dropping as (1 + z) 3, . However, as all the models at a given redshift 
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Figure 12. Flux decrement probability distributions in logarithmic bins of 
D at z = 2.4, 3, and 4 for the fiducial gas distribution, with inherent tem- 
peratures (black), at low resolution (grey), and with imposed H4 T-p (red). 
Also shown are the imposed H4 gas distribution with inherent temperatures 
(pink), with the imposed fiducial T-p (mustard), as well as both the fiducial 
and H4 gas distributions with an imposed flat T = 2 X 10 4 K (light blue 
and bright blue, respectively). The points are included to show the p(D) in 
the final (most opaque) bin. 



are normalized to the same (D), the effects of T-p changes on the 
mean absorption are removed. Though for visual clarity we do not 
plot the PDFs from the 2 x 1AA A particle simulation, the low reso- 
lution and fiducial PDFs differ by ~l-5%; this difference is much 
smaller than the typical model differences in Figures[TTland[T2l 
Thermal broadening reduces the number of transparent (D « 
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0) pixels, while pressure support increases the number of pixels 
at extreme flux decrements (for visual clarity, we do not plot the 
dark matter and non-fhermally broadened PDFs). At all redshifts, 
the spectra generated from the full H4 simulation have more pix- 
els with mid-range flux decrements, i.e. relative to the fiducial 
simulation, the H4 simulation produces fewer opaque pixels and 
fewer transparent pixels. This is as expected, since higher pressure 
smooths the gas distribution and higher temperature leads to greater 
thermal broadening. However, the fiducial simulation with the im- 
posed H4 T-p has fewer saturated pixels than the H4 simulation 
itself, and more pixels with D « 0.6-0.8. This surprising result 
highlights the sometimes complicated interplay between the den- 
sity distribution and the T-p relation. The fiducial simulation has 
more high overdensity gas (Figure |6), but because this gas is at 
high temperature (with T oc [1 + <5] ' 6 ), thermal broadening con- 
verts narrow, fully saturated features to broader, moderately satu- 
rated ones. 

The slope of the T-p relation has a direct impact on the flux 
PDF, independent of thermal broadening, because it affects the 
mapping from density contrast to flux. Equation l[8j gives r m oc 
(1 + S) 2 -°- 7a , so 7m oc (1 + S) 1S for our fiducial and H4 relations 
(with a « 0.6) and r„, oc (1 + S) 2 for our flat, T = 2 x 10 4 K rela- 
tion (a = 0). Imposing a constant T substantially alters the PDFs 
of both the fiducial and H4 simulations (blue curves in Figure [T2l. 
with the greater sensitivity of r m to (1 + S) leading to more satu- 
rated pixels and fewer mid-range pixels. The two simulations have 
significantly different PDFs, which with constant thermal broad- 
ening must arise from effects of pressure on the gas density dis- 
tribution. In contrast to the power spectrum and flux correlation 
function, thermal broadening and pressure support have compara- 
ble impact on the flux PDF, and they interact in complex ways. 
This different behavior arises because the PDF responds directly to 
the full overdensity distribution and its mapping to flux, while the 
power spectrum and correlation function measure the variance of 
this distribution as a function of scale. 



5 CONCLUSIONS 

We have investigated the relative importance of pressure support 
and thermal broadening in determining the longitudinal structure 
of the Lyman-a forest. Our main results come from comparing two 
SPH simulations with identical initial conditions but different pho- 
toionization heating rates, which produce a factor of ~ 2.3 dif- 
ference in the temperature of diffuse IGM gas. We have imposed 
different temperature-density relations on the simulation outputs to 
isolate physical effects, extracted spectra from the dark matter dis- 
tribution to extend our investigation to the pressureless case, and 
compared the fiducial simulation (2 x 288 3 particles) to a lower 
resolution simulation (2 x 144 3 ) to quantify numerical resolution 
effects. 

Equation l[3) shows that the Hubble flow across the Jeans 
scale is generally of the same magnitude as thermal broadening 
(HAj ~ oth). However, the IGM is an expanding, inhomogeneous 
medium evolving in the potential of a non-linear dark matter dis- 
tribution, so the Jeans length is at best an approximate description 
of the scale imposed by gas pressure support. It is therefore diffi- 
cult to know without detailed simulations whether thermal broad- 
ening or pressure support will dominate the structure of the for- 
est. A related question is the meaning of the density contrast S in 
the fluctuating Gunn-Peterson approximation (equation [8} ■ In prin- 
ciple, this should be the density of the gas that is absorbing Lya 



photons, but analyses of early SPH simulations found that Lya for- 
est spectra created from the (pressureless) dark matter distribution 
were remarkably similar to t hose created from the (pressure sup- 
ported) gas distribution (e.g. JCroft et alj[l998h . However, the rel- 
atively low resolution of those simulations (a mass resolution that 
is ~ 60 times lower than our fiducial simulation here) raised the 
possibility that both sets of spectra were artificially broadened to 
the numerical resolution limit. We investigate this issue more con- 
fidently here because the initial particle spacing of our 288 3 simu- 
lations, 43 hT 1 kpc comoving, is far below the typical Jeans length 
at IGM overdensities, Aj ~ 800/i _1 kpc comoving, and because 
our 144 3 simulation allows a direct resolution test. 

In broad brush, our conclusions are that thermal broadening 
dominates over pressure support in determining the visual appear- 
ance and statistical properties of the Lya forest, but that differences 
in gas pressure do have a noticeable effect. The widths of absorp- 
tion features are typically set by Hubble flow across the absorbing 
structure, though thermal broadening does smooth out small scale 
corrugations to create coherent features (see Fig. |4j- Once we in- 
clude thermal broadening, the spectra created from dark matter dis- 
tributions are similar to those created from the gas, even though the 
effects of pressure support on the gas density field are readily dis- 
cernible. Thus, one can drastically change the Jeans scale without 
drastically changing the forest. However, dark matter spectra are 
visually and statistically distinguishable from gas spectra, more so 
than in the earlier generation of lower resolution SPH simulations. 

Turning to individual statistics, we find that thermal broad- 
ening sets the turnover scale of the one-dimensional flux power 
spectrum, with the fiducial gas, H4 gas, and dark matter distri- 
butions producing similar power spectra if one imposes the fidu- 
cial or H4 T-p relation on all three. Similar conclusions hold for 
the coherence scale of the flux decrement autocorrelation function. 
However, in both cases, the weak impact of pressure support partly 
owes to a cancellation effect that arises with a sloped T-p relation: 
a higher pressure distribution has more "Jeans broadening," but it 
has less thermal broadening because there is less high overdensity, 
high temperature gas. When we impose a constant IGM tempera- 
ture of T = 2 x 10 4 K, the differences among the three cases are 
more noticeable, though they are still small compared to the effects 
of thermal broadening. For the flux decrement probability distribu- 
tion function, thermal broadening and pressure support have effects 
of comparable magnitude, though thermal broadening is still some- 
what more important. 

Our 144 3 and 288 3 simulations yield similar results for all our 
statistics. The resolution effects are larger at z = 4 than at lower 
redshifts, where they have a small but noticeable impact on all the 
statistics. For most purposes, the resolution of our 144 3 simulations 
(in a 12.5 ft" 1 Mpc comoving volume) is adequate. 

Though it is a stronger player than pressure support in set- 
ting the scale of the longitudinal Lya forest, thermal broadening 
is an inherently one-dimensional phenomenon. Because pressure 
acts in three dimensions, we expect it to play the main role in set- 
ting the transverse coherence of the Lya forest across neighboring 
sightlines. Growing samples of binary quasars with separations of 
A8 < 10" now make it po ssible to probe the expected Jeans scale 
dHennawi et al J2006l . [2009t) . We show in Paper II that the degree of 
transverse coherence on these scales is indeed sensitive to gas pres- 
sure support and insensitive to thermal broadening. Observational 
studies of close quasar pairs can directly probe the scale of pres- 
sure support on the Lyman-a forest, providing new insights into the 
physical state and thermal history of the high-redshift intergalactic 
medium. 
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